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A CRITICAL  ANALYSIS 

OF  CLIMATOLOGICAL  WIND  DATA  USED  IN  THE  FORECAST 
OF  RADIOACTIVE  DEBRIS  CLOUD  MOVEMENT 


I.  INTRODUCTION 


It  is  of  considerable  importance  to  communication  system 
performance  (ELF,  VLF,  HF)  in  a nuclear  environment  to  be  able  to  predict 
debris  cloud  transport  in  the  mesosphere.  Zalesak  and  Coffey  (1975) 
have  shown  that  transport  can  profoundly  alter  the  location  of 
radioactive  debris  clouds.  The  subsequent  beta  decay  within  the  cloud 
results  in  long  lasting,  widespread  ionization  that  can  severely  affect 
the  performance  of  communication  systems.  Accurate  prediction  of 
system  performance  depends  critically  on  our  ability  to  describe  the 
movement  of  debris  patches  by  mesospheric  wind  systems. 

The  purpose  of  this  report  is  to  assess  the  present  method  by 
which  the  spread  of  radioactive  debris  clouds  in  the  upper  atmosphere 
is  determined  and  to  suggest  guidelines  for  future  research.  The 
current  technique  used  to  forecast  debris  cloud  advection  uses  a 
Lagrangian  computer  code  with  model  wind  fields  ( Zalesak  and  Coffey, 
(1975)*  The  wind  fields  are  given  by  Groves  (1969),  CIRA  (1972).  These 
wind  fields,  being  based  upon  many  years  of  observational  data, 
represent  the  best  statistics  available  at  the  time  the  study  was 
undertaken.  The  results  produced  by  Zalesak  and  Coffey  clearly  hinge 
upon  the  accuracy  of  the  wind  field  data.  It  is  thus  important  to 
examine  the  data  carefully  and  ask  if  these  wind  fields  actually  give 
an  accurate  representation  of  the  upper  atmospheric  motions  relevant  to 
determining  debris  cloud  advection.  Groves'  data  are  an  "average"  or 
climatological  representation  of  the  wind  field  from  which  some  spatial 
and  temporal  fluctuations  have  been  removed  by  the  averaging  process. 

If  these  fluctuations  are  quite  small,  then  the  averaged  data  may  be  used 
to  accurately  forecast  the  transport  of  trace  constituents.  On  the  other 
hand,  if  the  fluctuations  are  large,  the  actual  wind  field  may  rarely 
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resemble  the  climatological  wind  field  and  resulting  debris  cloud 
forecasts  based  upon  the  latter  will  be  reliable  only  in  a climatological 
sense. 

In  this  report  the  wind  fields  given  by  Groves  shall  be 
examined  with  the  following  criteria.  First,  we  can,  to  some  extent, 
quantitatively  assess  variability  within  the  wind  field  data  by  examining 
the  standard  deviation  of  the  climatological  average  published  by  Groves 
(CIRA,  1972). 

Second,  we  examine  the  consistency  of  the  data  with  theoretical 
models  of  upper  atmospheric  dynamics.  While  inconsistency  between 
empirical  models  and  theoretical  models  does  not  necessarily  imply 
unreliable  data,  consistency  allows  us  to  use  theoretical  models 
where  empirical  data  may  be  lacking  or  difficult  to  obtain.  From  this 
viewpoint  we  can  determine  if  the  wind  structure  of  the  upper  atmos- 
phere as  computed  from  theoretically  postulated  heat  and  momentum 
sources  bears  any  resemblance  to  the  empirical  wind  structure  given 
by  Groves.  Or  conversely,  we  can  compute  the  implied  heat  and 
momentum  sources  required  to  maintain  Groves'  wind  model  and  compare 
with  known  sources.  Both  aspects  of  this  problem  will  be  discussed. 

Third,  the  zonal  wind  model  of  Groves  is  tested  for  stability 
to  small  wave  perturbations.  Instability  probably  implies  the  presence 
of  large  scale  eddy  mixing  which  could  greatly  affect  the  transport 
of  radioactive  debris. 

Within  the  body  of  this  report  the  data  and  the  observed 
variability  of  the  data  are  discussed  in  Part  II.  The  implied  heat 
and  momentum  sources  derived  from  Groves'  data  and  the  theoretically 
predicted  heat  and  momentum  sources  are  compared  in  Part  III. 

Stability  computations  for  Groves'  wind  model  are  presented  in  Part  IV. 

In  Part  V we  conclude  that  Groves'  wind  fields  are  probably  inadequate 
for  debris  cloud  advection  forecast  purposes  and  suggest  that  theoretical 
prediction  models  currently  under  development  can  be  used  to  provide 
more  reliable  results. 


II.  UPPER  ATMOSPHERIC  WIND  DATA 


Wind  observations  above  30  km  and  below  150  km  are  principally 
obtained  through  rocket  based  techniques.  At  high  altitudes  a 
meteorological  rocket  releases  an  object  or  chaff  which  is  tracked 
by  radar  as  it  falls.  Lateral  motion  of  the  falling  object  then 
yields  horizontal  wind  data  and  the  local  density  of  the  atmosphere 
may  be  computed  by  observing  the  rate  of  fall.  Compared  with 
radiosonde  observations  used  below  30  km,  rocket  methods  are  very 
expensive  and  technologically  complex.  As  a result,  the  network  of 
rocket  launching  stations  is  quite  sparse  and  regular  observations 
are  taken  only  weekly.  Figure  1 shows  the  station  locations  for  the 
Meteorological  Rocket  Network  (MRN) . 


The  rocket  data  obtained  through  the  MRN  facilities  contains 
both  systematic  and  random  deviations  or  "errors."  Both  kinds  reduce 
the  usability  of  the  derived  climatological  wind  models  for  forecasting 
We  may  further  subdivide  the  deviations  into  those  due  to  measurements 
(e.g.,  faulty  radar  techniques)  and  those  due  to  the  phenomena  (e.g., 
small  scale  eddies  inadequately  resolved  by  the  MRN  grid) . Quiroz 
(1969)  discusses  systematic  and  random  measurement  deviations  at 
length  and  his  findings  will  not  be  reviewed  here.  Measurement  error 
of  the  systematic  type  is  assumed  to  be  negligible,  whereas  random 
error  associated  with  the  measurements  is  assumed  to  be  removed  by 
climatological  averaging. 


Probably  the  most  obvious  source  of  systematic  deviations 
associated  with  phenomena  in  the  upper  atmosphere  is  that  produced 
by  the  presence  of  tidal  winds.  Most  MRN  data  is  taken  at  local  noon. 
Thus,  if  regular  diurnal  and  semidiurnal  tidal  wind  components  are 
present,  they  will  be  interpreted  in  any  local  climatological  analysis 
as  a component  of  the  mean  wind.  Lindzen  (1967)  has  computed  the 
amplitude  of  the  tidal  winds  up  to  100  km  and  has  found  that  winds 
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associated  with  the  solar  semidiurnal  and  diurnal  tides  may  be  as 
large  as  100  ms  ^ at  100  km  in  the  zonal  direction.  Measurement  of 
tidal  winds  below  60  km  shows  general  agreement  with  Lindzen's 
computation  with  some  disagreement  evident  above  60  km  (Glass  and 
Spizzichino,  1974). 

Groves'  wind  model  presented  in  CIRA  (1972) (also  Groves,  1969) 
has  been  constructed  by  grouping  MRN  and  other  data  into  monthly  or 
bimonthly  sets.  The  data  within  a set  have  been  further  subdivided 
into  four  hour  time  groups  depending  on  the  local  time  the  data  were 
taken.  The  average  within  each  group  was  computed  as  well  as  the 
mean  deviation.  Above  60  km  Southern  Hemisphere  data  were  assumed  to 
be  equal  to  Northern  Hemisphere  data.  Final  wind  model  values  were 
computed  by  an  iterative  scheme  involving  the  average  of  the  mean 
deviations  and  the  average  of  the  group  averages  and  a weighting 
formula.  Using  an  average  of  the  group  averages  is  equivalent  in 
some  sense  to  a daily  average.  Provided  large  monthly  changes  in  the 
amplitude  and  phase  of  the  tidal  components  do  not  occur  and  data 
samples  are  present  within  each  group,  this  method  should  eliminate 
the  systematic  error  introduced  by  tides.  In  reality,  however,  many 
groups  lack  data  altogether  above  60  km  so  that  model  points  are  based 
upon  only  one  or  two  groups  (cf . CIRA,  1972) . We  may  conclude  then 
that  high  altitude  winds  presented  by  Groves  probably  contain  con- 
siderable bias  from  tidal  winds  superimposed  upon  the  zonal  and 
meridional  mean  winds. 

Essentially,  the  systematic  deviation  introduced  by  tidal 
components  results  from  inadequate  temporal  resolution  of  the  zonal 
mean  wind  components.  Inadequate  spatial  resolution  can  also 
introduce  systematic  deviations.  In  particular,  quasistationary 
planetary  scale  waves  as  well  as  tides  have  wind  components  which 
vary  very  slowly  over  horizontal  distances  on  the  order  of  5,000  to 
10,000  km.  From  Figure  1 it  is  apparent  that  MRN  stations  are 
principally  located  in  the  northern  part  of  the  Western  Hemisphere, 


and  thus  the  network  will  be  unable  to  resolve  wind  components  associated 
with  very  long  zonal  scales. 

A comparison  of  West  European  and  North  American  data  presented  in 
CIRA  (1972)  indicates  the  presence  of  these  long  spatial  waves.  For 
example,  in  January  the  mean  zonal  wind  velocity  over  North  America  is 
^20  ms  * at  50  km  at  55°N,  while  the  mean  zonal  wind  velocity  over 
Europe  is  ^80  ms  ^ . The  difference  is  presumably  attributable  to  the 
long  wave  wind  components.  The  difference  between  North  American  and 
European  mean  zonal  wind  velocities  below  60  km  is  largest  during  winter 
and  is  consistent  with  the  observed  strength  of  planetary  scale  waves 
below  30  km  (van  Loon,  et  al.,  1973).  Theoretical  calculations  of  the 
amplitude  of  stationary  planetary  waves  in  the  upper  atmosphere  indicate 
that  these  waves  may  have  sizable  amplitudes  up  to  the  mesopause  and 
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may  generate  zonal  wind  components  as  large  as  20  ms  - 30  ms  in 
the  mesosphere  and  lower  thermosphere  (Schoeberl,  1975). 

Small  scale  eddies  are  probably  also  present  in  the  upper 
atmosphere  generated  by  baroclinic  instability  near  the  stratopause. 

If  their  horizontal  length  scales  are  much  smaller  than  1000  km, 
then  the  spatial  distribution  of  the  MRN  network  w’ill  be  inadequate  to 
properly  resolve  them.  These  eddies  will  appear  as  random  fluctuations 
in  the  MRN  data.  For  the  purpose  of  predicting  the  location  of  debris 
clouds,  these  eddies  may  be  as  important  as  the  zonal  mean  flow.  No 
information  is  available  from  Groves'  models  on  their  possible  struc- 
ture or  amplitude. 

All  of  the  phenomena  discussed  above  contribute  to  the  standard 
deviation  of  the  MRN  data  as  error.  In  Figure  2 the  model  values  of 
the  nean  zonal  wind  in  January  and  July  above  60  km  given  by  CIRA  (1972) 
and  the  standard  deviation  of  the  observations  from  the  model  values 
are  given.  Two  important  features  are  apparent.  First,  it  is  clear 
from  the  large  number  of  missing  standard  deviations  how  limited  the 
data  base  actually  is.  Second,  we  note  that  the  standard  deviation 
is  often  larger  than  the  mean  value  indicating  that  climatological  state 
of  the  wind  field  (Groves'  model)  occurs  as  an  exception  rather  than  the 
rule . 


III.  HEAT  AND  MOMENTUM  SOURCES  IN  THE  UPPER  ATMOSPHERE 


The  mean  zonal  circulation  is  driven  by  external  heat  and 
momentum  sources.  In  some  cases,  these  sources  may  be  theoretically 
computed  and  a circulation  model  developed  to  compare  with  observations 
(Leovy,  1964;  Baker  and  Strobel,  1975  , 1975,).  Alternatively,  a 
wind  model  derived  from  data  can  be  used  to  calculate  the  implied  heat 
and  momentum  sources  which  may  then  be  compared  to  theory  (Ebel,  1974). 
We  shall  consider  the  consistency  of  computed  wind  models  and  implied 
heat  and  momentum  sources  with  their  observed  and  theoretical  counter- 


in this  section 


Leovy  (1964)  showed  that  the  westerly  stratospheric  jet  observed 
in  the  winter  hemisphere  and  the  easterly  jet  observed  in  the  summer 
hemisphere  arise  from  the  meridional  ozone  heating  gradient  in  the 
stratosphere.  Mean  zonal  wind  maximums  of  80  ms  ^ were  computed  by 
Leovy  associated  with  mean  meridional  wind  velocities  of  0.7  ms  ^ . 

The  mean  zonal  wind  velocities  fluctuate  in  magnitude  throughout  the 
winter  (Belmont,  et  al. , 1975),  but  80  ms  * is  relatively  good  agreement 
with  Groves'  (1969)  climatological  value  considering  many  of  the 
simplifications  used  by  Leovy.  However,  Groves'  mean  meridional 
velocities  are  an  order  of  magnitude  larger  than  those  computed  by 
Leovy  and  Baker  and  Strobel.  Furthermore,  their  meridional  winds  blow 
from  the  summer  pole  to  the  winter  pole,  while  Groves'  meridional  winds 
are  quite  variable  depending  upon  latitude  and  altitude. 


The  discrepancy  between  these  computations  and  Groves'  data 
may  be  due  to  several  factors.  First,  assuming  Groves  meridional  winds 
are  correct,  the  mean  zonal  winds  (which  result  from  Coriolis  torques 
acting  upon  northward  or  southward  moving  flow)  may  be  computed  by  the 
following  equation. 


where  £2  is  the  earth's  frequency  of  rotation  and  0 is  the  latitude. 

£ is  the  Rayleigh  friction  coefficient;  v is  the  zonally  averaged 

meridional  velocity  of  the  wind,  and  u is  the  zonally  averaged  zonal 

velocity.  is  unknown  but  has  been  estimated  to  be  'VLO  ^ sec  ^ 

^ -1  

(Leovy , 1964).  Using  10  ms  for  v,  which  is  the  order  of  magnitude  given 
by  Groves  (1969),  gives  u '1000  ms  \ which  is  inconsistent  with  the 
TT  values  also  given.  If  IT  and  v are  assumed  correct,  we  are  forced  to 
conclude  that  equation  (1)  does  not  describe  the  correct  relationship 
between  u and  v,  and  the  addition  of  a momentum  source  term,  M,  of 
unknown  value  to  the  righthand  side  of  Equation  (1)  is  required  to  form 
a consistent  equation  between  u and  v.  It  is  also  apparent  that  the 
magnitude  of  M must  be  quite  large.  The  presence  of  eddies  which  could 
contribute  to  M are  known  to  exist  in  winter  but  are  generally  absent 
in  summer  (Kriester,  1972).  However,  large  values  of  v are  also 
indicated  by  Groves  for  the  summer,  so  this  explanation  is  implausible. 


A more  complete  calculation  of  the  required  momentum  and  heat 
sources  needed  to  maintain  the  Groves  model  winds  in  the  mesosphere 
(70  - 100  km)  has  been  carried  out  by  Ebel  (1974) . The  strength  of  the 
heat  sources  is  shown  in  Figure  3 for  solstice  conditions.  A 
comparison  with  the  computed  heat  sources  from  Park  and  London  (1973)  , 
Figure  4,  indicates  that  the  value  of  the  heat  sources  required  to 
maintain  the  Groves  wind  field  is  roughly  an  order  of  magnitude  too 
large.  Thus,  in  agreement  with  our  above  arguments,  it  is  improbable 
that  Groves'  v values  are  consistent  with  the  TT  values  given. 


Equation  (1)  is  the  zonal  mean  momentum  equation.  If,  as 
suggested  in  Section  II,  the  MRN  data  is  biased  by  tidal  and  stationary 
planetary  waves,  then  Equation  (1)  should  be  written  as 


where  (J>  is  the  geopotential  and  the  subscript  g indicates  the  values 
given  by  Groves  (1969)  which  are  now  not  assumed  to  be  zonal  means. 
Both  tidal  and  long  wave  components  can  theoretically  produce 


meridional  wind  velocities  as  large  as  those  reported  in  Groves' 

model  (Lindzen,  1967;  Schoeberl,  1975)  and  in  all  probability  it  is 

these  components  that  are  reported  by  Groves.  The  larger  velocities 

permitted  by  Equation  (2)  arise  from  the  presence  of  a zonal  pressure 

gradient  force  on  the  righthand  side  and  the  second  term  on  the  left- 

hand  side  which  is  an  inertial  term.  Both  of  these  tenns  are  much 

larger  than  the  term  B^u^  for  planetary  scale  anJ  tidal  motions.  The 

value  of  v is  thus  not  coupled  to  u alone.  For  tidal  and  planetary 
8 8 

scale  waves,  both  observation  and  theory  indicate  that  u 'W  '^10-20  ms 
in  the  stratosphere.  These  values  of  v are  more  consistent  with  the 
values  of  v^  and  suggest  that  the  data  are  indeed  biased  by  planetary 
wave  and  tidal  components. 


r 
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IV.  STABILITY  OF  MEAN  ZONAL  FLOW 


While  it  is  nearly  impossible  to  quantitatively  estimate  the 
magnitude  of  the  bias  that  long  wave  components  and  smaller  scale  eddy 
components  have  introduced  into  Groves'  wind  fields,  we  can  gain  some 
estimate  through  a stability  analysis.  Our  argument  is  as  follows: 

If  the  mean  zonal  wind  field  is  stable  to  wave  perturbations,  then 
any  finite  amplitude  eddy  disturbances  can  be  assumed  to  arise  from 
boundary  (tropospheric)  forcing.  If  the  flow  field  is  unstable,  then 
finite  amplitude  disturbances  may  arise  spontaneously  from  infinitesimal, 
local  disturbances. 

It  has  been  shown  by  Charney  and  Drazin  (1961)  that  only  large 
planetary  scale  eddies  can  propagate  into  the  upper  atmosphere. 

Synoptic  scale  disturbances  observed  in  the  troposphere  will  remain 
trapped  below  the  stratosphere.  Dickinson  (1973)  and  Simmons  (1975) 
have  shown  that  the  long  wave  components  are  the  fastest  growing  modes 
for  unstable  flow  fields  similar  to  those  observed  in  the  upper  strato- 
sphere and  mesosphere.  A computation  of  the  stability  of  the  observed 
zonal  mean  flow  field  as  given  by  Groves  (CIRA,  1972)  may  thus  indicate 
where  large  amplitude  eddy  components  could  arise. 

Using  the  Charney-Stern  stability  criteria  (Charney  and  Stern, 
1962)  we  compute  numerically  the  stability  of  Groves'  mean  zonal  wind 
field.  In  the  stability  criterion  for  an  atmosphere  bounded  by  rigid 
walls,  a necessary  condition  for  instability  is  that  Q,  defined  as 

2 — — 2 

_ _ . 9 u)  . _ _ q . 2 Q z9  e'9u)  . 

Q = 2(0+0))  - — r-  +-  3 tan  0 irw  - sin  0 e -5 r—  •=—  (3) 
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where  z = £n(pQ/p),  w = u/a  cos  0,  and  p is  pressure,  does  not  change 
sign  within  the  bounded  region. 
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Figures  5 and  6 show  the  value  of  Q computed  numerically  for 
parts  of  the  CIRA  (1972)  and  CIRA  (1965)  model  atmospheres,  respectively. 
Also  shown  is  a plot  of  the  corresponding  flow  pattern.  It  is  evident 
that  these  model  atmospheres  have  unstable  regions,  particularly  near 
the  stratopause  and  near  the  pole  at  all  levels,  as  indicated  by  the 
negative  values  of  Q.  A term  by  term  examination  of  Equation  (3) 
indicates  that  the  sign  change  in  Q is  principally  a result  of  sign 
changes  in  the  last  term.  Instabilities  developing  from  this  flow 
pattern  would  thus  be  primarily  baroclinic. 

The  assumption  that  Groves'  wind  models  characterize  the  mean 
zonal  flow  field  is,  of  course,  introduced  in  this  analysis.  One 
example  where  such  an  assumption  is  certainly  incorrect  is  evident  in 
Figure  4 where  a patch  of  negative  Q appears  to  50°N  and  30  km  in  the 
wind  model  taken  from  North  American  data.  No  such  region  appears  in 
the  computations  based  upon  European  data.  The  source  of  the  sign 
change  in  Q is  the  appearance  of  a region  of  easterly  winds  in  the 
midst  of  a westerly  jet  in  the  North  American  data.  The  winds  in  this 
region  are  probably  strongly  biased  by  planetary  waves  as  discussed  in 
Section  II  and  not  zonal  means;  hence,  the  Charney-Stern  stability 
criteria  does  not  apply. 

The  existence  of  high  frequency  motions  at  the  stratopause 
have  been  observed  by  Leovy  and  Akerman  (1973),  and  these  motions  may 
be  due  to  unstable  wind  configurations  such  as  those  given  by  the 
CIRA  (1972,  1965)  model  atmospheres.  In  any  event,  the  fact  that 
Groves'  models  are  unstable  indicates  that  either  eddy  components  of 
the  wind  have  biased  the  data  to  such  an  extent  that  the  wind  profile 
appears  unstable;  or,  that  eddy  components  are  present  with  sizable 
amplitudes . 
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V.  SUMMARY  AND  CONCLUSION 


The  usability  of  Groves'  climatological  wind  fields  for  fore- 
casting debris  cloud  advection  has  been  assessed  from  three  different 
viewpoints.  We  have  briefly  discussed  the  data  and  suggested  possible 
biasing  of  the  wind  field  by  tidal,  planetary  wave,  and  small  scale  eddy 
motions.  An  examination  of  the  standard  deviation  of  the  data  indicates 
that  the  actual  structure  of  the  wind  field  in  the  upper  atmosphere  rarely 
resembles  the  climatological  data  given  by  Groves. 

We  have  also  examined  the  heat  and  momentum  sources  implied  by 
wind  fields.  The  values  of  v given  by  Groves  (1969)  are  an  order  of 
magnitude  too  large,  and  imply  heat  and  momentum  sources  much  larger  than 
expected  from  theory.  We  conclude  the  v values  actually  represent  merid- 
ional velocities  associated  with  planetary  scale  waves  and  tides. 

Finally,  we  note  that  Groves'  zonal  wind  fields  are  unstable, 
especially  near  the  polar  stratopause.  The  instability  could  imply  the 
existence  of  large  amplitude  eddy  components  in  the  wind  field  for  that 
region.  We  conclude  that  while  Groves'  wind  fields  probably  represent 
the  best  available  data  they  nevertheless  inadequately  represent  the 
structure  of  the  upper  atmosphere  for  the  purposes  of  forecasting  the 
advection  of  debris  clouds. 

An  alternative  to  the  use  of  Groves'  wind  fields  to  forecast 
movement  of  a debris  cloud  is  the  use  of  a theoretical  prediction 
model  which  can  be  initialized  on  a regular  schedule  or  at  the  moment 
of  debris  cloud  release.  Such  a model  is  presently  used  in  the  lower 
atmosphere  and  gives  reliable  forecasts  up  to  three  days  in  advance. 

With  some  adaptations,  this  type  of  model  can  be  constructed  for  the 
upper  atmosphere  and  can  be  initialized  with  satellite  radiance  data. 

This  data,  which  is  in  the  form  of  temperature  fields,  is  currently 
available  up  to  60  km  (Chapman,  et  al.,  1972)  and  will  soon  be  available 
up  to  80  km  and  higher.  Upper  atmosphere  forecast  models  are  under 
development  at  NRL  at  present.  Madala.  et  al.,  (1975)  have  shown  that 
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tidal  winds  can  be  adequately  simulated  with  a spectral  forecast  model, 
and  Schoeberl  (1976)  has  been  able  to  determine  theoretically  the 
structure  of  planetary  waves  in  the  upper  atmosphere  using  a similar 
method. 
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Fig.  2 — January  and  July  statistical  wind  data  used  to  construct  the  CIRA  (1972) 
model  atmosphere.  Bold  face  numbers  show  the  standard  deviation  in  meters  per 
second  of  the  observational  data  at  the  indicated  altitudes.  X’s  indicate  only  a 
single  measurement  available;  — ’s  indicate  no  measurements.  Contours  labeled 
with  light  numbers  are  model  wind  values  in  meters  per  second. 


models  during  solstice  as  computed  by  Ebel  (1974).  Units  are 
0.864K  per  day.  The  winter  pole  lies  in  the  northern  hemisphere. 
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Fig.  5 — The  zonal  wind  model  (left)  and  value  of  Q defined  in  text  (right) 
computed  for  C1RA  (1972)  Northern  American  data  (top),  CIRA  (1972) 
West  European  data  (center),  and  all  longitudes  for  the  mesosphere  (bottom). 
Wind  profiles  for  January  were  used  in  all  instances.  Sign  changes  in  Q are 
necessary  conditions  for  instability  in  the  wind  profile. 
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Boeing  Company,  The 
P.  0.  Box  3707 
Seattle,  WA  98 124 

ATTN:  D.  Murray 

ATTN:  Glen  Keister 


B rown  Engineering  Company , The 
Cummings  Research  Park 
Huntsville,  AL  35807 

ATTN:  David  Lambert  MS  18 


California  at  San  Diego,  Univ  of 
Building  500  Mathews  Campus 
3172  Miramar  Road 
La  Jolla,  CA  92037 

ATTN:  Henry  G.  Booker 


Calspan  Corporation 
P.  0.  Box  235 
Buffalo,  N.Y.  14221 

ATTN:  Romea  A.  Deliberis 


Computer  Sciences  Corporation 
P.  0.  Box  530 

6565  Arlington  Blvd. 
Falls  Church,  V A 22046 
ATTN:  H.  Blank 

ATTN:  Barbara  F.  Adams 


Constat  Laboratories 
P.  0.  Box  115 
Clarksburg,  MD  20734 
ATTN:  R.  R.  Taur 


Cornell  University 

Department  of  Electrical  Engineering 
Ithaca,  N.  Y.  14850 

ATTN:  D.  T.  Farley,  Jr. 


EGSG,  INC. 

Los  Alamos  Division 

P.  0.  Box  809 

Los  Alamos,  N.M.  85544 

ATTN:  James  R.  Breedlove 


ESL , Inc. 

A95  Java  Drive 
Sunnyvale,  CA  93102 

ATTN:  J.  Roberts 

ATTN:  V.  L.  Mower 

ATTN:  James  Marshall 

ATTN:  R.  K.  Stevens 

General  Electric  Company 
Tempo-Center  for  Advanced  Studies 
8 1 6 State  Street  (P.  0.  Drawer  QQ) 

Santa  Barbara,  CA  93102 
ATTN:  Tom  Barrett 

ATTN:  Don  Chandler 

ATTN : DAS  I AC 

ATTN:  Warren  S.  Knapp 

General  Electric  Company 
Space  Division 
Valley  Forge  Space  Center 
P.  0.  Box  8555 
Philadelphia,  PA  19101 

ATTN:  M.  H.  Bortner,  Space  Sci.  Lab. 

General  Electric  Company 
P.  0.  Box  1122 
Syracuse,  N.Y.  13201 

ATTN:  F.  A.  Reibert 

General  Research  Corporation 
P.  0.  Box  3587 
Santa  Barbara,  CA  93105 
ATTN:  John  Ise,  Jr. 

Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 

ATTN:  Technical  Library 

ATTN:  Neil  Brown 

ATTN:  T.  N.  Davis 

GTE  Sylvania,  Inc. 

189  B.  Street 
ATTN:  Marshal  Cross 
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HRB-Singer,  Inc. 

Science  Park,  Science  Park  Road 

P.  0.  Box  60 

State  College,  PA  1 6801 

ATTN:  Larry  Feathers 

Illinois,  University  of 
Department  of  Electrical  Engineering 
Urbana,  IL  6 1 80 1 

ATTN:  K.  C.  Yeh 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 


ATTN: 

Ernest  Bauer 

ATTN: 

Hans 

Wol f hard 

ATTN: 

J.  M 

Aein 

ATTN: 

Joel 

Bengston 

Inti  Tel  & Telegraph  Corporation 
500  Washington,  Ave. 

Nut  ley,  N.  J.  07110 

ATTN:  Technical  Library 

ITT  Electro-Physics  Laboratories,  Inc. 

9140  Old  Annapolis  Road 
Columbia,  M0  21043 

ATTN:  John  M.  Kelso 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
8621  Georgia  Avenue 
Silver  Spring,  MD  20910 

ATTN:  Document  Librarian 

Lockheed  Missiles  6 Space  Co.,  Inc. 

P.  0.  Box  504 
Sunnyvale,  CA  94088 

ATTN:  Dept  60-12 

Lockheed  Missiles  S Space  Company 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

ATTN:  Billy  M.  McCormac,  Dept  52-14 

ATTN:  Martin  Walt,  Dept  52-10 

ATTN:  Richard  G.  Johnson,  Dept.  52-12 
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MIT  Lincoln  Laboratory 
P.  0.  Box  73 
Lexington,  MA  02173 

ATTN:  Mr.  Walden,  XI 13 

ATTN:  D.  Clark 

ATTN:  James  H.  Pannel 1 L-246 

ATTN:  Lib  A-082  for  David  M.  Towle 

Martin  Marietta  Corporation 
Denver  Distribution 
P.  0.  Box  179 
Denver,  CO  80201 

ATTN:  Special  Projects  Program  248 

Maxwell  Laboratories,  Inc. 

9244  Balboa  Avenue 
San  Diego,  CA  92123 

ATTN:  A.  J.  Shannon 

ATTN:  V.  Fargo 

ATTN:  A.  N.  Rostocker 

McDonnel  Douglas  Corporation 
530)  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
ATTN:  J.  Moule 

ATTN:  N.  Harris 

Mission  Research  Corporation 
735  State  Street 
Santa  3arbara,  CA  93101 
ATTN:  R.  Hendrick 

ATTN:  Conrad  L.  Longmire 

ATTN:  R.  E.  Rosenthal 

ATTN:  R.  Bougusch 

ATTN:  David  Sowle 

ATTN:  M.  Scheibe 

ATTN:  D.  Sappenfield 

ATTN:  P.  Fischer 

ATTN:  Ralph  Ki lb 
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Mitre  Corporation,  The 
Route  62  and  Middlesex  Turnpike 
P.  0.  Box  208 
Bedford,  MA  01730 

ATTN:  S.  A.  Morin  M/S 

ATTN:  Chief  Scientist  W.  Sen 

ATTN:  G.  Harding 

Mitre  Corporation,  The 
Westgate  Research  Park 
1820  Do  1 ley  Madison  Blvd. 

McLean,  VA  22101 

ATTN:  Allen  Schneider 

North  Carolina  State  Univ  at  Raleigh 
North  Carolina  State  Univ  Campus 
Raleigh,  N.C.  27507 

ATTN:  SEC  Officer  for  Walter  A. 

Pacific-Sierra  Research  Corp. 

1456  Cloverfield  Blvd. 

Santa  Monica,  CA  90404 

ATTN:  E.  C.  Field,  Jr. 

Philco-Ford  Corporation 
Western  Development  Laboratories  Div 
3939  Fabian  Way 
Palo  Alto,  CA  94303 

ATTN:  J.  T.  Matting  ley  MS  X22 

Photometries,  Inc. 

442  Marrett  Road 
Lexington,  MA  02173 

ATTN:  Irving  L.  Kofsky 

Physical  Dynamics,  Inc. 

P.  0.  Box  1069 
Berkeley,  CA  94701 

ATTN:  Joseph  B.  Workman 

RsD  Associates 
P.  0.  Box  3580 
Santa  Monica,  CA  90403 

ATTN:  Robert  E.  Lelevier 

ATTN:  Forest  Gilmore 

ATTN:  Richard  Latter 

ATTN:  William  B.  Wright,  Jr. 


I 


Flood 
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Rand  Corporation,  The 
1700  Main  Street 
Santa  Monica,  CA  90406 
ATTN:  Cullen  Crain 

Science  Applications,  Inc. 

P.  0.  Box  2351 
La  Jolla,  CA  92038 

ATTN:  Daniel  A.  Hamlin 

ATTN:  D.  Sachs 

ATTN:  E.  A.  Straker 

Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

ATTN:  L.  L.  Cobb 

ATTN:  Walter  G.  Chestnut 

A1TN:  Thomas  A.  Croft 

ATTN:  Donald  Neilson 

ATTN:  David  A.  Johnson 

ATTN:  Charles  L.  Rino 

ATTN:  E.  J.  Fremouw 

Stanford  Research  Institute 
306  Wynn  Drive, N.  W. 
Huntsville,  AL  35805 

ATTN:  Dale  H.  Davis 


Tri-Corn,  Inc. 

12215  Parklawn  Place 
Rockville,  MD  20852 

ATTN:  Darrell  Murray 


TRW  Systems  Group 
One  Space  Park 
Redondo  Beach,  CA  90278 
ATTN:  P.  H.  Katsos 

ATTN:  J.  W.  Lowry 

Visidyne,  Inc. 

19  Third  Avenue 

North  West  Industrial  Park 

Burlington,  MA  10803 


